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Abstract 

We present an elementary and explicit proof of the separability criterion 
for continuous variable two-party Gaussian systems. Our proof is based on an 
elementary formulation of uncertainty relations and an explicit determination 
of squeezing parameters for which the P-representation condition saturates 
the Sp(2, R) (g> Sp(2, R) invariant separability condition. We thus give the 
explicit formulas of squeezing parameters, which establish the equivalence of 
the separability condition with the P-representation condition, in terms of 
the parameters of the standard form of the correlation matrix. Our proof is 
compared to the past proofs, and it is pointed out that the original proof of the 
P-representation by Duan, Giedke, Cirac and Zoller(DGCZ) is incomplete. A 
way to complete their proof is then shown. It is noted that both of the 
corrected proof of DGCZ and the proof of R. Simon are closely related to our 
explicit construction despite their quite different appearances. 



1 Introduction 

The entanglement pQ is an intriguing property of quantum mechanics, but a quan- 
titative criterion of entanglement appears to be missing except for simple systems 
such as a two-spin system [21 E]. In view of this situation, it is remarkable that 
a quantitative sufficient condition for continuous variable two-party systems exists 
and that the criterion is necessary and sufficient for Gaussian states [U E]. The 
proofs given by Duan, Giedke, Cirac and Zoller (DGCZ) pE] and Simon [5] which 
consist of a series of logical steps are ingenious. However, their proofs are based 
on some specific notions and ideas in quantum optics, and thus their proofs are not 
readily accessible to those physicists who are interested only in the general aspects 
of entanglement in quantum mechanics. Moreover, these two proofs are seemingly 
quite different and their mutual connections are not clear. This problem and related 
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issues have been discussed in the past by several authors [6j [7J [8j [9j [TOj [TTJ [TJl [13] • 
The present status of the quantum separability problem is nicely reviewed in [2]. 

We here present an elementary and explicit proof by starting with the elemen- 
tary analysis of Heisenberg uncertainty relations in the manner of Kennard [T5| [T6] 
and an explicit determination of squeezing parameters which establish that the P- 
representation condition saturates the Sp(2, R) <g> Sp(2, R) invariant separability 
condition. We thus give the explicit formulas of squeezing parameters, which estab- 
lish the equivalence of the separability condition and the P-representation condition, 
in terms of the parameters of the standard form of the correlation matrix (or second 
moments). It is also pointed out that the original proof of the P-representation by 
DGCZ is incomplete, and a way to complete their proof is shown. It is then shown 
that both of the corrected proof of DGCZ and the seemingly quite different proof 
of Simon are closely related to our explicit construction. 

Our treatment is based on a clear recognition that the separability condition 
associated with uncertainty relations is invariant under general Sp(2, R) <g> Sp(2, R) 
transformations, whereas the condition for the P-representation of Gaussian states 
is not invariant under general Sp(2, R) <g> Sp(2, R) transformations. A combination 
of these two apparently contradicting relations is the basis of our construction of 
the explicit solution. 

2 Entanglement and Kennard's relation 
2.1 Kennard's relation 

We consider a two-party system (or a two-particle system in one- dimensional space) 
described by canonical variables (qi,pi) and (q 2 , Pi)- We define 

X(d, f) = diqi + d 2 pi + hq 2 + /2P2, 

X(g, h) = g x q x + g 2 pi + hq 2 + h 2 p 2 (2.1) 

where all the coefficients 

d T =(d 1 ,d 2 ), f T =(f u f 2 ), g T =(g u g 2 ), h T = (h u h 2 ) (2.2) 

are real numbers. The Kennard's relation for a mixed state p = Y2k -Pfcl^fcXV^I with 
Pk > and Y2k — 1 is written as [S] (for any choice of d ~ h) 

((AX(d, f)) 2 ) + ((AX(g, h)f) > \d T Jg + f T Jh\ (2.3) 

where we defined 

((AX(dJ)) 2 )=Tr{(AX(dJ)) 2 p} (2.4) 
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with AX(d,f) = X(d,f) - (X(d,f)) and (X(d,f)) = Tr{X(d,f)p}, for example, 
and the 2x2 matrix 

The relation (2.3) is derived from Tr{f}tfp} > and Trjr^p} > for fj = 
AX(d, f)+iAX(g, h), and the right-hand side of (2.3) stands for | [AX(d, f), AX(g, h)]\. 
We examine (2.3) more precisely by starting with 

<(AX(rf,/)) 2 > = ((x(dJ)-(X(d.f))) 2 ) 

= ^P fc ((x(rf./)-(X(rf./))) 2 ) fe 

k 

= p k({x(d.f) - (X(d, f)) k + (X(d, f)) k - (X(d, f)}) 2 } k 
k 

= J2 p k[((x(d.f) - (x(d, f)) k ) 2 ) k + (<x(d, f)) k - (x(d, /))) 2 ] 

k 

> ^P fc ((x(rf,/)-(X(rf,/)) fe ) 2 ) fc (2.6) 

k 

which holds for a general mixed state for any real numbers d and /. The equality 
sign holds only for 

(X(dJ)} k -(X(dJ)}=0 (2.7) 
for all k where (X(d, /)) = E fe P k (X(d, f)) k with 

(X(d, f)) k = J d gi dq 2 r k (qu Q2)X(d, f)M<li, fc). (2.8) 

The condition (2.7) is trivial for a pure state, but it imposes a stringent condition 
on a mixed state. The Kennard's relation for a general pure state is given by (2.3) 
if one sets P k — 1 for specific k and others zero 

{(X(d,f) - {X{d,f)) k ) 2 ) k + ((X(g,h) - (X(g,h)) k ) 2 ) k 
> \d T Jg + f T Jh\ (2.9) 
for any d ~ h, which is also written as 

t 2 {(x(d,f) - {X(d,f)) k ) 2 ) k + {(x(g,h) - (X(g,h)) k ) 2 ) k 

-t\d T Jg + f T Jh\>0 (2.10) 



by replacing d td, f — > tf for any real t and thus the discriminant gives the 
conventional form of Kennard's relation. The Kennard's relations for pure states 
imply 

J2Pk[((x(d,f) - {X(d,f)) k ) 2 ) k + ((x(g,h) - {X(g,h)) k ) 2 ) k } 



k 

>\d T Jg + f T Jh\ (2.11) 

for any d ~ h, which is more precise than (2.3) because of the removal of extra terms 
as in (2.6). 

2.2 Separability condition 

For a separable pure state ip k {qi,q2) = (pkiq^fk^), we have 

((x(d,f)-(X(d,f)) k y) k = ((d l q 1 + d 2 p l -(d 1 q 1 + d 2 p l ) k ) 2 ) k 

+ <(/!& + hn - (/ig 2 + hV2)k)\ (2.12) 
and similarly for (^X(g, h) — (X(g, h))kj ) k - We thus have 

[((x(d,f) ~ {X(d,f)) k ) 2 ) k + ((x(g,K) - (X{g,h)) k ) 2 ) k ] 
= [{(X(d,0) - (X(d,0)) k y) k + ((X{g,0) - (X(g,0)) k y) k 

+((x(o,f) - (x(o,/)) fc )V + {(x(o,h) - (x(o,h)) k y) k \ 

> \d T Jg\ + \f T Jh\ (2.13) 

which holds for any d ~ h. Here we used (2.9) for / = h — or d — g — 0. The 
equality sign holds only for 

[(diqi + d 2 pi - (diqi + d 2 pi) k ) + % (giqi + g 2 pi - (#i<7i + g 2 Pi) k)}<Pk(qi) = 0, 
+ /2P2 ~ (fifa + /2P2>fc) + i (hfa + h 2 p 2 - (hiq 2 + h 2 p 2 ) k )](f k (q 2 ) = 0, 

(2.14) 

for suitable d ~ h with d T Jg > and f T Jh > 0. 

We finally conclude from (2.6) and (2.13) for any separable density matrix 

((AX(rf,/)) 2 ) + ((AX(^/ i )) 2 ) 

> Pk[((X(d, f)) k - (X(d, J)))' + ((X(g, h)) k - (X(g, hrf] 
k 

+ \d T Jg\ + \f T Jh\ (2.15) 



for any d ~ h. 

We next define the variables (£ a ) = (qi, p±, q 2 , p 2 ) and the 4x4 correlation matrix 
Vby 

V=(V af) ), V af 3 — ^{A£ a A£p + A£pA£ a ) — ^({A| a , A^}) (2.16) 
with A£ a = £ a — (£ a ), which can be written in the form 

v = ( b ) < 2 - 17 > 

where A and -B are 2x2 real symmetric matrices and C is a 2 x 2 real matrix. We 
also define 

V = (V aP ), V a/3 = J2 P k(Ma)k(Mfi)k (2.18) 

k 

and 

V- = ( £ i ) (2^9) 

where A and B are 2x2 real symmetric matrices and C is a 2 x 2 real matrix. 
Both of V and V are non-negative. This quantity V plays a central role in the 
P- representation. 

The basic relation (2.15) for separable states is then written as 

d T Ad + fBf + 2d T Cf + g T Ag + h T Bh + 2# T C/i 

> d T Ad + fBf + 2d T Cf + # T i# + /i T 5/i + 2g T Ch 

+ \d T Jg\ + \f T Jh\. (2.20) 

while the Kennard relation for general states is written as 

d T Ad + fBf + 2d T Cf + g T Ag + /i T 5/i + 2g T Ch 

> d T Ad + fBf + 2d T C/ + £ T i# + h T Bh + 2# T C7i 

+ |rf T J^ + / T J/i|. (2.21) 

Note the difference between \d T Jg\ + \f T Jh\ and |<i T Jg + f Jh\. 
The antisymmetric commutator parts in 

(AlaAfc) = i({A| a , A&}> + i([A| a , A&]) (2.22) 



which may be added to A and B in (2.17) do not contribute to (2.20) since d ~ h 
are all real. 

Under the S^^G Sp(2, R) ® Sp(2, R) transformations of (<?i,Pi) and (^2,^2), 
respectively, we have 

A^SxAgf, B^S 2 BS%, C^S x CSl (2.23) 

which is equivalent to the transformation 

d-Sfd, /->S 2 T /, /^S 2 T /i (2.24) 

in (2.20) if one recalls J = S\JSf, J = S^JS^ ; the inequality (2.20), which 
is valid for any d ~ /i, holds after the transformation (2.24) and in this sense 
(2.20) is invariant under the above Sp(2,R) ® Sp(2,R). To be precise, we do not 
use any property of the wave function under Sp(2,R) ® Sp(2,R), and thus our 
5p(2, i?) (g) 5p(2, i?) transformation is rather defined by (2.23) for given constant 
matrices A, B and C. 

The difference between the two expressions in (2.20) and (2.21) appears when 
one replaces B and C by SjBS 3 and CS 3 (and also B and C by SjBSs and CS 3 ), 
respectively, with 



^3=(J \ (2.25) 

One can undo the replacements in the first expression (2.20) by transformations 
/ ~^ £3/ an d h S^h, whereas it leads to {dfjg — f T Jh\ in the second Kennard 
relation (2.21). The separability condition thus demands that the Kennard relation 
should hold both for the original system and for the system with the replacements 
of B and C by SJBS3 and CS3, respectively, which may a priori be unphysical for 
inseparable systems. By using S 3 , one can adjust the signature of detC at one's 
will 0. 

It is also useful to consider the separability condition (2.20) with subsidiary 
conditions g = J T d and h = J T f, 

d T Ad + f T Bf + 2d T Cf + d T JAJ T d + f T JBJ T f + 2d T JCJ T f 

> d T Ad + f T Bf + 2d T Cf + d T JAJ T d + f T JBJ T f + 2d T JCJ T f 

+ {d T d + ff) (2.26) 

and with subsidiary conditions g = J T d and h = —J T f 

d T Ad + f T Bf + 2d T Cf + d T JAJ T d + f T JBJ T f - 2d T JCJ T f 

> d T Ad + f T Bf + 2d T Cf + d T JAJ T d + f T JBJ T f - 2d T JCJ T f 

+ (d T d + f T f). (2.27) 



6 



For general inseparable states in (2.21), we have the first condition (2.26) only if 
one wants to keep (d T d + f T f) on the right-hand side in this form. The basic 
Sp(2, R) <g> Sp(2, R) invariance of uncertainty relations (2.20) and (2.21) is lost in 
these conditions (2.26) and (2.27) with subsidiary conditions, but they have appli- 
cations in the analysis of the P-representation. 



3 Separability and P-representation in Gaussian 
states 

3.1 General analysis 

One can bring any given V in (2.17) by Sp(2, R) <S> Sp(2, R) transformations to the 
standard form [U [5] (see also Appendix A) 



V 
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Cl 


o \ 
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c 2 
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b 





\ 





C-2 





b J 



(3.1) 



One may understand the relations (2.20) and (2.21) ( and also (2.26) and (2.27)) 
in two different ways: 

(i) In the first interpretation, one may understand these relations (2.20) and (2.21) 
as an infinite set of uncertainty relations (and their variants) for any given constants 
d ~ h. In this interpretation, the relations (2.26) and (2.27) correspond to the ones 
used by DGCZ jl] if one chooses d and / suitably. 

(ii) In the second interpretation of the relations (2.20) and (2.21), one may under- 
stand these relations holding for any choice of d ~ h and thus imposing constraints 
on the allowed ranges of the elements a, b, ci, c 2 of the standard form of Vq in (3.1), 
for example. We adopt this second interpretation, which was also adopted by Si- 
mon [5]. In this interpretation, the full relation (2.20) is more restrictive than the 
relations (2.26) and (2.27) with the subsidiary conditions on d ~ h. In other words, 
the elements a,b,c±,C2 which satisfy (2.20) automatically satisfy (2.26) and (2.27), 
but not the other way around. In our analysis below, we interpret these relations as 
constraints on Ci, c 2 for fixed a, b. 

The separability criterion is given by (2.20). On the other hand, the P-representation 
depends on the condition (see Appendix B) 

V > \l (3.2) 
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namely 



d T Ad + f T Bf + 2d T Cf > X -{d T d + ff) (3.3) 



for any d ~ /. By using a special property of the P-representation, namely, a special 
property of the coherent state, one can also write (3.2) as 

P-=(£ i)=V-i/>0 (3.4) 

or 

d T Ad + f T Bf + 2d T Cf 

= d T Ad + / T 5/ + 2d T C/ - i (d T d + / T /) > (3.5) 

for any d ~ /. Here P _1 agrees with V in (2.19). See Appendix B. 

We first note that the P-representation implies the separability condition, since 
(3.5) also implies 

g T Ag + h T Bh + 2g T Ch 

= g T Ag + h T Bh + 2g T Ch - -{g T g + h T h) > (3.6) 

2 

and thus adding these two relations (3.5) and (3.6), we have 

d T Ad + f T Bf + 2d T Cf + g T Ag + h T Bh + 2g T Ch 
= d T Ad + fBf + 2d T Cf + g T Ag + h T Bh + 2g T Ch 

l -(d T d + f T f)+ l -{g T g + h T h). (3.7) 

When one combines this relation with 

l -(d T d + f T f)+ l -(g T g + h T h) > y/(<Fd + frf)(gTg + hTh) 

= ^{d T d + f T f)(g T J T Jg + h T J T Jh) 

> \d T Jg\ + \f T Jh\, (3.8) 

one reproduces the separability condition (2.20). This is natural since the P- 
representation is in fact separable. 

But the inverse is not obvious. The separability condition (2,20) is invariant 
under Sp(2, R) ®Sp(2, R) in (2.23), whereas the P-representation condition (3.2) or 
(3.3) is not invariant under Sp(2, R) ® Sp(2, R) by noting that 

d T SfS 1 d + f T S^S 2 f + d T d + ff (3.9) 
8 



in general. In this sense these two conditions cannot be equivalent to each other. 
One may write the condition for P-representation as 



V-^SS T >0 (3.10) 

for a suitable but arbitrary S G Sp(2,R) ® Sp(2,R); this relation implies the or- 
dinary P-representation condition S~ 1 V(S~ 1 ) T — \l > for a suitable Sp(2,R) ® 
Sp(2,R) transformed S~ 1 V(S~~ 1 ) T . Written in the form (3.10), the condition for 
P-representation has a formally invariant meaning in the following sense. For any 
St G Sp(2, R) ® Sp(2, R), we have 

SiVSf " \s 1 SS T Sf > (3.11) 

which is written as 

V - ^S'(S') T > (3.12) 
with V = StVSf amd S' = S ± S G Sp{2, R) ® ^(2, R). 



3.2 A new explicit proof 

We here present an explicit proof of the separability criterion for continuous variable 
two-party Gaussian systems. Our explicit construction gives the formulas of squeez- 
ing parameters, which establish the equivalence of the separability condition with 
the P-representation condition, in terms of the parameters of the standard form of 
the correlation matrix (3.1). 

When one regards the separability condition as a constraint on the range of |ci| 
and | C2 1 in the standard form V$ (3.1), it is written as 

A(ab - c\){ab - c 2 2 ) > (a 2 + b 2 ) + 2|c x c 2 | - i, 

v / (2a- 1)(26- 1) > |ci| + |c 2 | (3.13) 

together with a > 1/2 and b > 1/2. The conditions a > 1/2 and b > 1/2 are 
respectively derived by setting / = h = and d = g = in (2.20). The first relation 
in (3.13), which was derived by Simon [5], corresponds to 

4det[V + ^(o (3 - M) 
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up to a transformation S3 in (2.25), and thus it is manifestly invariant under 
Sp(2,R) <S> Sp(2,R). The second condition in (3.13) is derived by the weaker con- 
ditions in (2.26) and (2.27) for the standard representation Vq, and it is used to 
exclude the nonsensical solutions of (3.13) with c 2 — > 00 and c\ — > 00 for fixed a and 
b. The conditions (3.13) are equivalent to (2.20). 

The separability condition (3.13) is explicitly solved as 

cl < ^{[2a6(l + t 2 ) +t]- 2^/a 2 b 2 (\ - t 2 ) 2 + t{a + bt){at + &)}, 



c\ < i{[2a6(l + t 2 ) +t}- 2^a 2 b 2 {\ - t 2 ) 2 + t(a + bt)(at + b)} 



(3.15) 



for 



0<t= \c2\l\d\ < 1 



(3.16) 



where we choose | C2 1 < \ci\ without loss of generality. 
On the other hand, one may choose S in (3.10) as 



S(r 1 ,r 2 )S T (r 1 ,r 2 ) = 



( i/n \ 

On 

l/r 2 

V r 2 J 



(3.17) 



with suitably chosen r± > 1 and r 2 > 1. By choosing the standard form of Vq in 
(3.1), the eigenvalues of Vq — \S '(ri, r 2 ) S* 7 (ri, r 2 ) are given by 



< A >>± = - + " - ± v ((a - - {b ~ k )? + 4c?} - 

(A 2 ) ± = j{(o - ir,) + (6 - ir 2 ) ± ^((a - \r,) - (b- \r 2 )Y + 4cl}(3.18) 

The P-representation exists if (Ai)± > and (A 2 ) ± > 0, namely, if the following two 
conditions are simultaneously satisfied 



1 1 

(a - -n)(6- -r 2 ) > cl 



(3.19) 



together with a>\, b > \, (a - ±) + (b - ^) > 0, and (a - |n) + (6 - ±r 2 ) > 0. 



2r 2 - 
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When one regards (3.15) and (3.19) as constraints on the pair of variables 

(cl,c 2 2 ) (3.20) 

for given a and b, the P-representation condition is more restrictive than the sepa- 
rability conditions, namely, the set of points {c\, c|) allowed by the P-representation 
condition (3.19) always satisfy the separability condition (3.15). To be precise, we 
are working on the line defined by t 2 = c 2 /c\. We thus expect that these two con- 
ditions can coincide only for the extremal value of the P-representation condition 
(3.19) with respect to r\ and r 2 with fixed t. We show that this is indeed the case. 
We thus want to prove 

(«-^)(* 



2r/ v 2r 2 (t,r 1 ) / 
= ^[(a-^)(&-^r 2 (f,ri))] (3.21) 

= ^{[2a6(l + t 2 ) +t\- 2^/aW(l - t 2 ) 2 + t(a + bt)(at + b)} 

for a suitable 1 < r\ < 2a (and 1 < r 2 < 2b) for any given < t < 1 by regarding r 2 
as a function of r\ and t. By this way we establish that the separability condition 
(3.15) agrees with the P-representation condition (3.19) with a suitable Sp(2, R) <g> 
Sp(2, R) transformation. 

We start with the equality in the left-hand side of (3.21) 

and take the derivative of the both hand sides with respect to r x with fixed t. We 
then have 

\—(b - 1 ) + ( a - 1 ) 1 dr2 ] 



which is solved as 

dr 2 _ (fe-W^^ + ^-s^fey) 
dr i ( a -^i) + ^iy( a -^) 

We next consider the stationary point (or extremal) of 



(3.24) 



(a-^i)(6-ir 2 (t,n)) (3.25) 
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with fixed t, namely 

1 /, 1 . . . 1 , 1 <9r 2 

--^-- 2 r^r 1 ))-(a--r l) -- r ^ (3.26) 

This relation combined with (3.24) gives rise to 

(« - - = (a - \r l)( i - (3.27) 

The relations (3.22) and (3.27) give 



n = ^T7i, r 2 = f^- (3.28) 



which are symmetric in r\ and r 2 . The relations in (3.28) are solved as 

n = ——{ab(l - t 2 ) + ^a 2 b 2 (l - t 2 ) 2 + t(a + bt)(at + &)}, 
at + b 

r 2 = — *— - {ab(l - t 2 ) + Ja 2 b 2 {l - t 2 ) 2 + t{a + bt)(at + b)} (3.29) 
a + bt 

with 0<t=|c 2 |/|ci|<l, which determine the squeezing parameters. 
We see from (3.29) that 

n = r 2 = 1 (3.30) 

for t — 1, and 

n = 2a, r 2 = 26 (3.31) 

for t = 0. One can also confirm 

ri a + at 

oo > — = > 1, 

* -a6(l - t 2 ) + v/a 2 6 2 (l - t 2 ) 2 + t(a + at) (at + 6) ~ 

oo > — = , > 1 3.32 

t -ab(l - t 2 ) + v /a 2 6 2 (l - t 2 ) 2 + t(a + bt)(at + 6) ~ 

by noting t(at + o) < (a + at) and t(a + at) < (at + 6) for < t < 1 and the triangle 
inequality. By recalling (3.28), we thus conclude 

2a > n > 1, 26 > r 2 > 1 (3.33) 
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for a > | and b > ~. 

We finally evaluate by using ri and r 2 in (3.29) 



1 . ,a + 6t, Ja 2 lP{l - t 2 ) 2 + t(a + bt)(at + b). 

— — [a + at{ -) — 

At 2[ at + b J at + b 



rj , ,at + b, Ja 2 b 2 (l - t 2 ) 2 + t(a + bt)(at + b). 
x[b + bt{- 1 



' a + bt a + bt 

1 



^{[2afe(l + t 2 ) +t]- 2^/a 2 b 2 (l - t 2 ) 2 + t(a + bt)(at + b)} (3.34) 

which is a remarkable identity. This relation establishes (3.21), namely, the fact 
that the boundaries of the conditions for separability and P-representation coincide 
for any < t = | C2 1 / [ Ci | < 1. 

Our explicit construction proves that the P-representation condition (3.19) with 
suitably chosen S'(r 1 , r 2 ) G Sp(2, R) <g> Sp(2, R), where 1 < r x < 2a and 1 < r 2 < 2b, 
is equivalent to the separability condition (3.15) for any < t = |c 2 |/|ci| < 1, 
and thus the separability condition (3.15) is a necessary and sufficient separability 
criterion for two-party Gaussian systems. Our formulas of r\ and r 2 in (3.29) give 
the explicit expressions of squeezing parameters to achieve the above equivalence in 
terms of the parameters of the standard form of the correlation matrix (3.1). 

4 Comparison with the past proofs 
4.1 Proof of Duan,Giedke,Cirac and Zoller 

The analysis of Duan, Giedke, Cirac and Zoller (DGCZ) [1] starts with the constraint 



'22-1) (--1) 



(2ar l - 1) (26r 2 - 1) 

which is written as 



(4.1) 



^-1) (--1) 

(4.2) 



(rari — 1) (mr 2 — 1) 

by noting 2a = n and 2b = m in their notation of the Sp(2, R)®Sp(2, R) transformed 
correlation matrix 

/ nri c-s/TxT-i \ 

n/r 1 d / ^r 1 r 2 



M 



c^/t\T2 mr 2 

\ d I y/rir 2 m/r 2 J 



(4.3) 
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with n > m > 1 and |c| > |c'| > 0. Their normalization corresponds to M = 2V. 

One can rewrite (4.2) as r 2 (mr 2 — 1) = X(ri){m — r 2 ) and solve this quadratic 
equation in r 2 in the form 



, , 1 - X ± J(l - X) 2 + Am 2 X , , 

r 2 (n)± = 2 m (4.4) 



where we defined 



X(n) = |L=i> = rs^l) , 4 . 5) 

which assumes X(l) = 1, X(n — e) = oo, X(n + e) = — oo and X(oo) = — oo. Here 
e is an infinitesimal positive quantity which is eventually set to 0. One thus finds 

r 2 (l) + = l, r 2 (l)_ = -l, 

r 2 (n — e) + = m — e, r 2 (n — e)_ = — oo, 

r 2 (n + e) + = oo, r 2 (n + e)_ = m + e, 

r 2 (cx)) + = oo, r 2 (oo)_ = m (4.6) 

and thus the solution has a rather involved branch structure 0. 



1 By writing (4.5) as 

X n ee i = - n(n - ») + — - (2n 2 - 1) 

n — ri ri — n 

one can show — oo < X(ri) < — [n + \/n? — l] 2 for n < ri < oo and the upper bound is achieved 
at r\ = n + V^ 2 — 1. Since X > for 1 < n < n, one can confirm that the content inside the 
square root in (4.4) 

(1 - Xf + Am 2 X = (X + (m + \J m 2 - l) 2 ) (x + (m - \Jm 2 - l) 2 ) 

is positive definite for 1 < n < oo when n > m > 1. Consequently the solutions r2{r\)± in (4.4) 
are real, and r^iri)-^ for 1 < n < n and r2(ri)- for n < ri < oo defines a continuous real function 
?"2(ri) for 1 < ri < oo, which is all that is necessary in the analysis in [3]. Incidentally, the content 
inside the square root above becomes negative and thus r 2 (ri)± become complex for m > n for 
some interval in n < r% < oo. In particular for n — m, r 2 — T\ is a solution. 
When one rewrites (4.4) as 

2m 



l-i± v /(l — T^^T 



for X > by taking X inside the square root and continues this expression for negative X also, 
then the single branch Ti(r\).\. covers the entire domain 1 < n < oo for n > m. This is because 
one picks up an extra — sign when one takes negative X inside the square root. 
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One may next consider [4] 



-[V(nrx-l)(mr 2 -l)-./(- - 1)(- - 1) ] (4.7) 

V T\ r 2 

which satisfies (if one chooses the first branch r 2 (r 1 ) + in (4.4)) 

/(I) = \c\ - |c'| > 0. (4.8) 

They then use the bound [4J 



1 / m + 1 

c < \ n(m ) = \ n(m — 1) , (4-9) 

V m V 

which is derived from (2.26) and (2.27) by setting g = 0. One may now examine 
in (4.7) in the domain n < r\ < oo [I] by choosing the branch r 2 (ri)_ in (4.4) 
for which one can show m < r 2 (ri)_. One then establishes /(oo) < by using (4.9) 
as in [1]. 

One thus concludes a solution for /(ri) = in the interval 1 < r\ < oo, as shown 
in [4], for the set of states which satisfy the condition (4.9). Incidentally, 



f(n) = y/nm\c\ jLL - ^J{n 2 - l)(m 2 - 1), (4.10) 

and thus if one can establish 

^fnm\c\ ^ — — < a/ (n 2 — l)(m 2 — 1), (4-H) 



|c 



one then has /(n) < 0. So far we briefly summarized the analysis in [4] together 
with a comment on (4.10) and (4.11) which are used later. 

We now examine the main issue of the separability condition, eq.(16) in [4], 

2 ni + n 2 mi + m 2 2 1 

fl ° 2 2^ |ci| - |c 2 | > a + ^ (4-12) 



with a 2 , = w ^3^" — \/ ^TT' w hich gives rise to 



vVx - l)( Wl - 1) ± v/(m 2 - l)(n 2 - 1) > |ci| + |c 2 | (4.13) 
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for 1 < ri < n and n < r 1; respectively, if one recalls that n 2 < 1 and m 2 < 1 for 
n < ri when (4.3) is written in the form 



M 



( 


m 







o \ 







n 2 





c 2 




Cl 










\ 





c 2 





m 2 ) 



(4.14) 



Note the appearance of the crucial ± sign in (4.13) because of (n 2 — 1) ' 1712 1 



ri2 — 1 



n 2 



^1 V rta-l = ~~ ( m2 — -'-)( n2 — -0 f° r n 2 < 1, for example. It appears that 
this minus sign was overlooked in [4]. 

If /(ri) = has a solution in the interval 1 < r\ < n, one can derive the condition 
for the P-representation (eq.(17) in [4]) 

VK-llK-l) > | Cl |, v/(m 2 -l)(n 2 -l) > |c 2 | (4.15) 

by combining the first relation in (4.13) with /(ri) = 0, and the proof of the P- 
representation in [4] naturally goes through. On the other hand, if f{r{) = has a 
solution in the interval n < r\ < oo, one finds that the separability condition (4.13) 
with n < r± is inconsistent with /(ri) = for |c 2 | 7^ since one then has 



- 1)K - 1) - v 7 ^ - l)(n 2 - 1) = \ci\ - |c 2 | > |c x | + |c 2 |. (4.16) 

This puzzling result for < r x may indicate that some essential information is 
missing to analyze the P-representation. This is indeed the case as shown below. 
The condition M — / > of the P-representation in fact requires 

(m - 1) + (mi - 1) > 0, (n 2 - 1) + (mj - 1) > (4.17) 

in addition to (4.16), since the eigenvalues of M — I are given by 

(Ai)± = \[{n x - 1) + (m x - 1) ± sj ((m - 1) - (mi - l)) 2 + 4c?], 

(A a )± = l[(n 2 - 1) + (m 2 - 1) ± ^/((n 2 - 1) - (m 2 - l)) 2 + 4c|]. (4.18) 

See (3.18) in a different representation. It is obvious that if (4.17) are not satisfied, at 
least one of the eigenvalues in (4.18) becomes negative. More intuitively, M — I > 
cannot be maintained for sufficiently small c\ and c 2 if (4.17) are not satisfied. One 
can confirm that (4.15) is derived from the requirement (Ai j2 )± > in (4.18), which 
is equivalent to M — / > 0, only under the conditions (4.17). These extra conditions 
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(4.17) are missing in [4], and the condition (ji2 — 1) + (m 2 — 1) > is violated in the 
case with n < r\ < oo for which n 2 < 1 and m 2 < 1, and thus no P-representation 
exists for n < r% < oo 0. 

It has been argued in [4] that separability or inseparability is independent of 
squeezing, but when one replaces separability by the P-representation one finds that 
the P-representation is very sensitive to squeezing. No information about (4.17) is 
contained either in the separability condition (4.12)(eq.(16) in [4] in terms of EPR- 
like operators) or in their analysis of /(ri) = 0, and in this sense one may conclude 
that the proof the P-representation in the original scheme of [4] is incomplete. 

One may now recall that the exact separability condition (3.14) is Sp(2,R) Cg> 
Sp(2, R) invariant while the P-representation condition is not invariant as is shown 
in (3.9). The separability condition (3.15) is thus independent of squeezing param- 
eters while the P-representation condition (3.19) explicitly depends on squeezing 
parameters. The squeezing is an auxiliary device to show that the P-representation 
condition combined with suitable squeezing is equivalent to the Sp(2, R) <g> Sp(2, R) 
invariant separability condition. A salient feature of the analysis in |4j is that they 
use the separability condition (4.12) which depends on squeezing parameters. 

To prove the P-representation starting with the separability condition, one needs 
to satisfy two conditions (4.17) and (4.15). A way to achieve this purpose in the 
framework of [4] is to show (4.15) by using (4.12) for the squeezing parameter in the 
range 

1 < n < n. (4.19) 

We then automatically ensure (4.17) by means of (4.2) (and (4.6)) which is always 
assumed in our analysis. The non-trivial task is to show (4.15). For this purpose, we 
start with the first relation in (4.13), which is derived from the separability condition 
(4.12) if (4.2) is satisfied for 1 < r\ < n. It is then confirmed that the first relation 
in (4.13) with t\ = n (and thus r 2 = m because of (4.2)) 

I c > I 

Jim 2 - l)(n 2 - 1) > Icil + |c 2 | = s/rvm\c\ + -±=L (4.20) 

y'nm 

ensures f(n) < in (4.10). We thus conclude that there exists a solution for 
f{ r i) — in the interval 1 < r\ < n for the state which satisfies the separability 

2 It appears that for most cases in practice, the continuous variable states automatically satisfy 
Lemma 2 (standard form II) in [4], namely the matrix M in (4.3) with the constraint (4.2) for a 
suitable n in the interval 1 < r\ < oo which satisfies f(r\) — in (4.7). It should be emphasized 
that standard form II in 0] as it stands is based on the quite weak condition (4.9) and thus it is valid 
for inseparable states also. The standard form II holds even for n < r*i but the P-representation 
does not exist for n < r\ and m < r2, namely, for ri2 < 1 and m,2 < 1. 
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condition (4.12) for any r\ in 1 < r\ < n. The proof of the condition M — I > for 
the P-representation is then complete. Namely, (4.15) (eq.(17) in |4j) together with 
(4.17) is established by combining the first inequality in (4.13) with f(rx) = for 
any |c| > \c'\ in (4.3). An important new ingredient of the present scheme compared 
to the original scheme in [4] is that the order of the analyses of /(ri) = and the 
separability condition (4.12) is reversed and the separability condition (4.12) with 
1 < ri < n now plays a central role in the analysis of f(r\) = 0. One can confirm 
that (4.15) is equivalent to (3.19) when converted into our notation. 

We here add further comments on the scheme of DGCZ in view of our explicit 
construction in Section 3. 

Firstly, it is interesting that their condition (4.2), of which origin is not clearly 
stated in [4], agrees with our extremal condition (3.27). 

Secondly, it is shown that the weaker forms of the separability condition, (2.26) 
and (2.27), when applied to the representation (4.14) give rise to the condition 



V[( n i + n 2) ~ 2] [(mi + m 2 ) 
To be explicit, (2.26) gives 



2] > |ci| + |c 2 |. 



(4.21) 



/ n 1 + n 2 

n l +n 2 

ci + c 2 

\ ci + c 2 



Cl + c 2 


mi + m 2 




\ 

Cl + c 2 


mi + m 2 / 



21 > 



by taking M = 2V into account, and (2.27) gives 



/ n\ + n 2 




V 



ci - c 2 




ci - c 2 
ni + n 2 

mi + m 2 
-Ci + c 2 



\ 

-ci + c 2 


mi + m 2 J 



21 > 0. 



One can also show 

v /[{n 1 + n 2 )-2][(m l +m 2 )-2] > y/[ ni - l][mi - 1] + ^ [n 2 - \}[m 2 - 1} (4.22) 

where the equality holds only when the condition (4.2) is satisfied. This relation 
(4.22) is established by considering 



f(x) = a/Ni + x(n 2 - m)) - l][mi + x(m 2 - mi)) - 1] 



(4.23) 
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with 

f"(x) = A[( ni -l)(m 2 -l)-(n 2 -l)( mi -l)} 2 

x [mi + x(m 2 - mi)) - Vf' i/2 [ni + x(n 2 - ni)) - 1]~ 3/2 < 

(4.24) 

except for (4.2), namely, 

{n 2 - 1) = (m 2 - 1) 
(m - 1) (mi - 1) 

for which f'(x) = 0. By using the property of the convex function 2/(1/2) > 
/(l) + /(0) one can establish (4.22); the condition 1 < ri < n is sufficient to keep 
f(x) real for < x < 1. Their separability condition (4.13), which is derived from 
(4.21) when the equality in (4.22) holds, thus corresponds to the weaker form of the 
separability condition. 

This fact suggests that under the extremal condition (4.2), the weaker forms 
of the separability condition, (2.16) and (2.27), are sufficient to ensure the P- 
representation if supplemented by an additional constraint 1 < r± < n. 



4.2 Proof of Simon 

The analysis of the case c\c 2 > by Simon [5] is quite elegant. Starting with the 
standard form of Vq in (3.1) and applying a set of Sp(2, R) Cg> Sp(2, R) transforma- 
tions, he arrives at the form of V 



V 



( ay 2 x 2 


ciy 2 




c lV 2 \ 

a/(y 2 x 2 ) c 2 /y 2 

by 2 /x 2 

c 2 /y 2 bx 2 /y 2 ) 



which is also written as 



V 



ar l c Xs Jr\r 2 \ 

a/n c 2 j^r x r 2 

C\\/r\r 2 br 2 

c 2 j^r^F 2 b/r 2 ) 



by defining 



(4.25) 



(4.26) 



n = (xy) 2 , r 2 = y 2 /x 2 



(4.27) 
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He uses the crucial condition 

Cl c 2 



ax 2 — b/x 2 a/x 2 — bx 2 



(4.28) 



which allows the diagonalization of V by a Sp(4, i?) transformation. This Sp(4, R) 
preserves the Kennard relation, and thus one can use the Kennard relation to show 
the P-representation. 

The condition (4.28) is written as 



i r 1 cia + c 2 b a + (c 2 /ci)b 



r 2 c 2 a + Cib {c 2 /Ci)a + b 



which implies (for a > b) 



T\ a 

1 < — < r- 4.30 

r 2 6 



It is interesting that the condition (4.29) agrees with our explicit construction (3.29). 

Simon [5] shows that V in (4.25) can be diagonalized by an Sp(4, R) transfor- 
mation as V = diagf _) with 



iy 2 {ax 2 + 6/x 2 ± ^(ax 2 - 6/^ 2 ) 2 + 4c 2 }, 
-?T 2 {a/x 2 + 6x 2 ± y^a/x 2 - bx 2 ) 2 + Ac 2 }. (4.31) 



2 

The equality of two smaller eigenvalues = is ensured if 

4 a/x 2 + bx 2 — Jialx 2 — bx 2 ) 2 + 4c 2 

y = r x r 2 = 

ax 2 + b/x 2 — y (ax 2 — b/x 2 ) 2 + 4c 2 



i + b(r 1 /r 2 ) - Z( a ~ K r i/ r Z) 2 + 4c|(ri/r 2 ) 
i(ri/r 2 ) + 6 - y/{a(r 1 /r 2 ) - b) 2 + 4c 2 (ri/r 2 ) 



(4.32) 



and the Kennard's relation (uncertainty relation), K-k'_ > 1/4, implies k_ = k/_ > 
1/2 which in turn gives rise to 

( ar i ~ \)iP r 2 ~ ^) > c 2 rir 2 , 

(a/ri^)(^^)>4/(nr 2 ) (4.33) 
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or equivalently 



lr\ 2r 2 



(a - -n)(b - -r 2 ) > c\. 



(4.34) 



This last relation naturally agrees with the condition of the P-representation (3.19). 
The separability condition, which agrees with the Kennard's relation for cic 2 > 0, 
thus ensures the P-representation. The case cic 2 < is equally treated by replacing 
Ci and c 2 by |ci| and |c 2 |, respectively [5]. 

Although we have not succeeded in proving (4.32) by using our explicit solution 
(3.29) due to technical complications, one can confirm that (4.29) combined with 
one of the relations (3.28) gives rise to our explicit solution (3.29). We thus believe 
that the solution given by Simon agrees with our explicit construction. 

5 Conclusion 

We have presented an elementary and explicit analysis of the separability criterion 
of continuous variable two-party Gaussian systems. In particular, we derived the 
explicit formulas of squeezing parameters, which establish the equivalence of the 
separability condition with the P-representation condition, in terms of the param- 
eters of the standard form of the correlation matrix (or second moments). In the 
course of our analysis, we corrected the shortcomings of the past proof of DGCZ [I]. 
Our explicit construction also clarified the basic equivalence of the past seemingly 
quite different proofs of the separability criterion [4J [5] in the sense that both of the 
past proofs are closely related to the present explicitly constructed solution. 

I thank K. Shiokawa for an informative discussion. 

A Standard form of V 

We recall the elements of Sp(2, R) 




(A.l) 
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which satisfy SJS T = J. One can bring V in (2.17) to the standard form 



V 



( 


a 





Cl 


o \ 







a 





c 2 




Cl 





b 





\ 





c 2 





b > 



(A.2) 



by suitable Sp(2, R) <g> Sp{2, R) transformations [HE]; real symmetric A and B can 
be made diagonal by two-dimensional rotations with suitable parameters 9 in (A.l) 
and then applying the second elements in (A.l) with suitable parameters x, A and 
B are made proportional to the unit matrix. After these transformations C remains 
real. By applying a suitable two-dimensional orthogonal transformation S\ ® S2, 
which is an element of Sp(2, R) ® Sp(2, R), we can diagonalize C 



D\Ls D2 



ci 
c 2 



(A.3) 



By this way we arrive at (A.2). 



B P-representation 

We define the generating function of all the correlations (or moments) of dynamical 
variables by 

X(X, rj) = Tr(pexp{i(Aigi + \ 2 pi + 771^2 + »)}) (B.l) 

where Ai ~ r]2 are real numbers. By expanding xi^iV) i n powers of Ai ~ r]2, one 
can generate all the moments of dynamical variables. Following the convention in 
this field, we define the Gaussian states by 

X(A, r]) = exp{-i(Ai, A 2 , Vi,V2)V(\ 1 , A 2 , 771, r? 2 ) T } (B.2) 

where V is the correlation matrix in (2.17), namely, all the correlation functions are 
determined by the second moments. One can write (B.l) as 

X (A, n) = Tr(pexp{i(A*a + Aa f + rj*b + rjb 1 )}) (B.3) 

with 

a = 4=(gi + ipi), b= -^=(q 2 + ip 2 ), 

A = ^|(Ai + «A 2 ), v = -^=(ri 1 + iij 2 ) (B.4) 
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The Gaussian state is called P-representable if the density matrix is written as 

p = J d 2 a J d 2 l3P(a,[3)\a,[3)(a,[3\ (B.5) 

where \a, (3) is the coherent state defined by 

a\a, (3) — a\a, (3), b\a, (3) — (3\a, (3), (a, (3\a, (3) — 1 (B.6) 
or to be explicit 

\a,(3) = e a£l '-^ 2 \0}®eP~ bt -^ 2 \0}. (B.7) 

Thus the P-representable states are separable. 

By using the density matrix (B.5) in (B.3) and after normal ordering the expo- 
nential factor in (B.3), we have 

X(\v) = f d 2 a J d 2 (3P(a,[3)exp{i(\*a + \a* + r]*f3 + r]f3*)} 

xexp{-i(|A| 2 +M 2 )} (B.8) 
or, if one combines this expression with (B.2) we have 

ex P{— ^0i> A 2 , 171,172)0" ~ ^ J )(^i^2,^i,^) T } 

= J d 2 a J d 2 (3P(a, (3) exp{i(Ai«i + A 2 a 2 + r^fa + V2P2)} (B.9) 

with a = (a 1 + ia 2 ) /V2 and p = (p ± + ifo) /V2. Thus P(a, (3) in (B.9) is given by 

\fdetP 1 

P(a,(3) = ^^e W {--(a 1 ,a 2 ,(3 1 ,(3 2 )P(a 1 ,a 2 ,(3 1 ,(3 2 ) T } (B.10) 

where 

P~ 1 = V --I>0 (B.ll) 

2 - v I 

which defines the condition for the P-representation. 

The formula (B.9) indicates that the right-hand side generates the correlations 
of the form 

J d 2 a J d 2 (3P(a, [3)(a, (3\tf\a, [3)(a, [3\a\a, [3), (B.12) 

for example, which may be compared to (2.18). By recalling (B.5), this relation 
shows that all the second moments in the right-hand side of (B.9) are given by V 
in (2.18) (if one chooses (a 1 ") = (a) = (¥) = (b) = 0). This property establishes the 
special relation (3.4) of the P-representation. 
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